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Abstract 

This paper presents some formulae to calculate moments of inertia for solids of revolution and for 
solids generated by contour plots. For this, the symmetry properties and the generating functions of the 
figures are utilized. The combined use of generating functions and symmetry properties greatly simplifies 
many practical calculations. In addition, the explicit use of generating functions gives the possibility of 
writing the moment of inertia as a functional, which in turn allows to utilize the calculus of variations 
to obtain a new insight about some properties of this fundamental quantity. In particular, minimization 
of moments of inertia under certain restrictions is possible by using variational methods. 

PACS {45.40.-F, 46.05.th, 02.30.Wd} 
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1 Introduction 

The moment of inertia (MI) plays a fundamental role in the mechanics of the rigid rotator £Q, and is a 
useful tool in applied physics and engineering Hence, its explicit calculation is of greatest interest. In 
the literature, there are alternative methods to facilitate the calculations of MI for beginners 0], at a more 
advanced level [5], or using an experimental approach [Hj. Most textbooks in mechanics, engineering and 
calculus show some methods to calculate Mi's for certain types of figures Q][2|[3|. However, they usually do 
not exploit the symmetry properties of the object to make the calculation easier. 

In this paper, we show some quite general formulae in which the Mi's are written in terms of generating 
functions. On one hand, by combining considerations of symmetry with the generating functions approach, 
we are able to calculate the Mi's for many solids much easier than using common techniques. On the other 
hand, the explicit use of generating functions permits to express the Mi's as functionals; thus, we can use 
the methods of the calculus of variations (CV) to study the mathematical properties of the Ml for several 
kind of figures. In particular, minimization processes of Mi's under certain restrictions arc developed by 
applying variational methods. 

The paper is organized as follows: in Sees. [2HH]we derive expressions for the Mi's of solids of revolution. 
In Sec. 0]we calculate Mi's of solids by using the contour plots of the figure, finding formulae for thin plates 
as a special case. Sec. |3J shows some applications of our formulae, exploring some properties of the Ml by 
using methods of the CV. Sec. El contains the analysis and conclusions. 

2 Mi's for solids of revolution generated around the X-axis 
2.1 Moment of inertia with respect to the X-axis 

Let us evaluate the moment of inertia of a solid of revolution with respect to the axis that generates it, 
in this case the X— axis according to Fig. ^ We shall assume henceforth, that the solid of revolution is 
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Figure 1: Solid of revolution generated from the X—axis. The unit of integration is a piece of hoop with 
cross sectional area (dx) (dr x ) and length r x dO. We show a rotated view of the piece of hoop in the upper 
right corner, for an easier visualization of its dimensions. 

generated by two functions f\ (x) and fi (x) that fulfill the condition < f\ (x) < $2 (x) for all x £ [xq, £/]■ 
Owing to the cylindrical symmetry exhibited by the solids of revolution, it is more convenient to work 
in a cylindrical system of coordinates. The coordinates are denoted by x,r x ,9 where r x is the distance from 
the X—axis to the point, the coordinate 9 is defined such that 9 = when the vector radius is parallel to 
the Y — axis, and increases when going from Y (positive) to Z (positive) as shown in Fig. ^ 

We consider a thin hoop with rectangular cross sectional area equal to (dx) (dr x ) and perimeter 2nr x as 
Fig. ^displays. Our infinitesimal element of volume will be a very short piece of this hoop, lying between 9 
and 9 + dO, the arc length for this short section of the hoop is r x dO. Therefore, the infinitesimal element of 
volume and the corresponding differential of mass, are given by 



the distance from the X—axis to this element of volume is r x . Therefore, the differential moment of inertia 
for this element reads 

dl x = r 2 x dm = rip (x, r x ,6) (dx) (dr x ) (dO) , 
and the total moment of inertia is given by 

r x f ( rh(x) \ i-e f ~] 



In this expression we first form the complete hoop (no necessarily closed), it means that we integrate in 9 
because when we form the hoop, the x and r x variables are maintained constant and the integration only 
involves 9 as a variable. After the completion of the hoop, we form a hollow cylinder of minimum radius 
fi (x), maximum radius fi (x) and height dx. We make it by integrating concentric hoops, where the radii 
r x of the hoops run from f\ (x) to fi (x). Clearly, the variable x is constant in this step of the integration. 
Finally, we integrate the hollow cylinders to obtain the solid of revolution, it is performed by running the x 
variable from xq to Xf as we can see in Fig^ This procedure gives us the formula of Eq. I|2.2[l . 

The expression given by Eq. (|2.2|) is valid for any solid of revolution (characterized by the generating 
functions f\ (x), fi (x)) which can be totally inhomogencous. Even, the revolution does not have to be from 
to 27T. If we assume a complete revolution for the solid with p — p(x, r x ) the formula simplifies to 



dV = (dx) (dr x ) r x (d9) ; dm = p (x, r x ,9) dV , 



(2.1) 




(2.2) 




(2.3) 
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and even simpler for p = p (x) 



p(x) f 2 (x) 4 -f 1 (x) 4 



dx 



(2.4) 



We see that all these expressions for the MI of solids of revolution along the axis of symmetry, are written 
in terms of the generating functions fi (x) and f 2 (x). In particular, it follows from Eq. (|2.4[) that for 
homogeneous solids of revolution (and even for inhomogeneous ones whose density depend only on x i.e. the 
height of the solid), the volume integral involving the calculation of MI is reduced to a simple integral in 
one variable. We point out that common textbooks misuse the cylindrical properties of these type of solids, 
evaluating explicitly all the three integrals even for homogeneous objects. 



2.2 Moments of inertia with respect to the Y and Z axes 

Y 




Figure 2: Distance from the Y—axis to a differential element of volume. In terms of the coordinates defined 
in Sec. \2.1\ the square distance yields r 2 = x 2 + r 2 sin 2 9. 

Let us calculate the moment of inertia of the solid of revolution with respect to the Y — axis. We use the 
same element of volume of the previous section. The square distance from the Y — axis to such element of 
volume is x 2 + r 2 sin 2 8, as shown in Fig. [21 Therefore, the moment of inertia of this element of volume with 
respect to Y reads 

dI Y = (r 2 sin 2 6 + x 2 ) dm , 
with dm given by Eq. 12. 1|) . Integrating in a way similar to the previous section, the MI becomes 



Iy = 







fSf 




f"\ 




/ p(x,r x ,9)sm 2 d6 


r 3 




Jh{x) 


Je 






p(x,r x ,6) d9 



f x dr x > x dx 



(2.5) 



Once again, this formula is valid for any inhomogeneous solid of revolution (even incomplete) generated by 
the functions f\ (x), fi (x). When assuming a complete solid of revolution with p = p(x,r x ), and taking 
into account Eq. 1)2.3(1 . the latter formula reduces to 



Iy = + 2n 



p(x,r x ) r x dr x 



dx 



(2.6) 



From this expression we derive the interesting property 7y > which is valid for any complete solid of 

revolution with azimuthal symmetry. Additionally, for p = p (x) we have 



Iy = Ij^Ix + 7T 



p(x) x 2 f 2 (x) - f\ (x)' 



dx 



(2.7) 
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By the same token, for the Z— axis, which is perpendicular to the X and Y— axes and with the origin as 
the intersection point, we have the following general formula 



Iz 











f"\ 




1 p(x,r x ,9)cos 2 9 dO 


r 3 

' X 


!x 


J hi*) 


Je 







rMx) 








f"\ 




/ p{x,r x ,9 


) dO 






Jfi(x) 


J9 








(2.8) 



We emphasize that in the case of a complete revolution, the formula (|2.8|) . coincides exactly with Iy in Eq. 
(|2.6|) when p docs not depend on 9, as expected from the cylindrical symmetry. Indeed, for Iz to be equal 
to Iy, the requirement of azimuthal symmetry could be softened by demanding the conditions 



p = p(x,r x )p(9) ; / p{9)sm z 9d9 
Je 

for if the conditions l|2.9() are held, we get that 



p (9) cos 2 9 d9 , 



(2.9) 



p (9) cos 2 9 d9 



p(9) d9 , 



00 



and Iy = Iz even for an incomplete solid of revolution with no azimuthal symmetry. 

From Eqs. (|2. 512. 8(1 . we see that for the calculation of the MI for axes perpendicular to the axis of 
symmetry, we use the same limits of integration as for the symmetry axis; thus, we do not have to care 
about the partitions. Once again, these Mi's are written in terms of the generating functions of the solid 
fi (x) and f 2 (x). 

Finally, we emphasize that textbooks do not usually report the moments of inertia for solids of revolution 
with respect to axes perpendicular to the axis of symmetry. However, they are important in many physical 
problems. For instance, a solid of revolution acting as a physical pendulum requires the calculation of such 
Mi's, see example |HJ 




X 



Figure 3: Frustum of a right circular cone with a conical well as a solid of revolution. The shadowed surface 
is the one that generates the solid. /i (x) and f 2 (x) are the functions that provide the limits of integration. 



Example 1 Mi's for a truncated cone with a conical well (see Fig. 0). The generating functions read 

AW = 

/a (a) = [^TT^jx + ai, (2-10) 



fl(l-f) if xe[0,h] 
if x e (h, H] 

a\ — a 2 ' 



II 
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where all the dimensions involved are displayed in Fig. For uniform density, we can replace Eqs. }2.1U\) 
into Eqs. {^J ^TJjjj to get 



'X 



TTf) 

To 



{H (af + a\ + a±al + a\a 2 + a\a 2 2 ) - R A h) 



Iy = Iz = -^Ix 



npH 3 



1 



-aia 2 + a x 



6 



R 2 
6 \H 



It is more usual to give the radius of gyration (RG) instead of the MI. For this we calculate the mass of the 
solid by using Eq. \A.1\) . finding 



irp 
T 



M='^&[H (oia 2 + a\ + a\) - R 2 h] 



(2.11) 



The radii of gyration become 



3 {H (af + a\ + a^l + a\a 2 + a\a 2 2 ) - R 4 h} 
10 [H (aia 2 + al+al) - R 2 h] 



Ktr = Ki = 



K 



x 



-H 



1 n n J- „2 , 1-2 R 2 ( /l N 



2 5 [H (oi0 2 + a? + al) - R 2 h) 



(2.12) 



By making R = (and/or h = 0) we find the RG's for the truncated cone. With R = and a\ = 0, we gei 
the RG's of a cone for which the axes Y and Z pass through its base. Making R = and a 2 = 0, we find 
the RG's of a cone but with the axes Y and Z passing through its vertex. Finally, by setting up R = 0. and 
o-i = a 2 ; we obtain the RG's for a cylinder. In many cases of interest, we need to calculate the Mi's for axes 
XcYcZc passing through the center of mass (CM), these Mi's can be calculated by finding the position of 
the CM with respect to the original coordinate axes, and using Steiner's theorem (also known as "the parallel 
axis theorem"). Applying Eqs. ^A.S^A.^]) the position of the CM for the truncated cone with a conical well 
is given by (xcm,0, 0) with 

[(2oiaa + 3a? + a 2 ) H 2 - R 2 h 2 } 



XCM = 



4 [H (aia 2 



R 2 h 



(2.13) 



Gathering Eqs. \2.12L \2.13jH we find 



K 



K 



x i 



K 2 



Ki 
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2.3 Another alternative of calculation and a proof of consistency 

In addition to the the parallel and the perpendicular axis theorems, there is another useful theorem about 
Mi's that is not usually included in common texts, namely \7\ 

I X +lY + Iz = ^Y. m ^ R2 ^ ( 2 - 14 ) 

i 

where X, Y, Z are three mutually perpendicular intersecting axes, m% is the mass of the i— th particle and R{ 
is its distance from the intersection. We shall see that our general formulae for Mi's of solids of revolution, 
fulfill the theorem. From Eqs. (JO ED wc nave 

Ix+Iy + Iz = 2 I ' ' ] f f (x 2 + r 2 ) p (x, r x ,9) r x (d9) {dr x ) (dx) . (2.15) 

Jx Jfi{x) Je 

Moreover, if we take into account that the distance from the intersecting point (the origin of coordinates) to 
the element of volume is R 2 = x 2 + r 2 , and using Eq. (|2.1fl we conclude that 

Ix + Iy + Iz = 2 / R 2 dm , (2.16) 
Jv 
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which is the continuous version of the theorem established in Eq. I|2.14[) . As well as providing a proof of 
consistency, this theorem could reduce the task to estimate the Mi's, especially when a certain spherical 
symmetry is involved. 

Further, it is interesting to see that the Mi's Ix,Iy,Iz, in Eqs. i|2.2l 12.51 12.8fl fulfill the triangular 
inequalities 

Ix < Iy + Iz , 

and same for any cyclic change of the labels. The triangular inequalities follow directly from the definition of 
MI, and are valid for any arbitrary object. Though the demostration of these inequalities is straightforward, 
they are not usually considered in the literature. In the case of thin plates, one of them becomes an equality. 

Example 2 The following example shows the usefulness of the theorem of Eqs. \2.1J\ \ 2.1 ftp in practical 
calculations. Let us consider the MI of a sphere centered at the origin, whose density is factorizable in 
spherical coordinates such that p — p(R). Where R is the distance from the origin of coordinates to the 
point. The symmetry of p leads to Ix — Iy = Iz and the theorem in Eq. $2.16)) gives 

pRo 

3I X = 8tt / p{R) R 4 dR , (2.17) 
Jo 

the mass of the sphere is 

,6 

M = Air / p{R)R 2 dR , (2.18) 
Jo 

from which the moment of inertia can be written as 



2 f n p(R)R 4 dR 
Ix = — A/ , ; . 2.19 

We can calculate for instance, the classical MI of an electron in a hydrogen-like atom, with respect to an axis 
that passes through its CM. For example, for the state (1,0,0) we have that p(R) = 2 (Z/a^ff^ 2 e ~ ZR / a ° j 
where Z is the atomic number and Oq the Bohr's radius. The MI becomes 

8m e 2 
l x - ~^T a o ■ 

3 Mi's of solids of revolution generated around the Y-axis 




Figure 4: Solid of revolution generated around the Y—axis. The unit of integration is a piece of hoop with 
cross sectional area (dy) (dr y ) and length r y dep. 

By using a couple of generating functions f\ (x) and fi (x) like in the previous section, we are able to 
generate another solid of revolution by rotating such functions around the y—axis instead of the X— axis 
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as Fig. 21 displays. In this case however, we should assume that x$ > 0; such that all points in the 
generating surface have always non- negative x coordinates. Instead, we might allow the functions f\ (a;), 
f 2 (x) to be negative though still demanding that fi (x) < f 2 [x) in the whole interval of x. In this case, it 
is more convenient to use another cylindrical system in which we define the coordinates <j)), where r y 

is the distance from the Y— axis to the point, and the angle </> has been defined such that = when the 
vector radius is parallel to the Z— axis (positive), increasing when going from Z (positive) to X (positive). 
One important comment is in order, since the surface that generates the solid lies on the XY— plane, the x 
coordinate of any point of this surface (which is always non-negative according to our assumptions) coincides 
with the coordinate r y , therefore we shall write fi (r y ) and f 2 (r y ) instead of fi (x) , f 2 (x) for the functions 
that bound the generating surface. 

The procedure to evaluate the MI in the general case is analogous to the techniques used in sectional 
the results arc 



i-x 





' rh{r v ) 


f"\ 

1 Xq 


Jh{r y ) 



4>} 




p(r y ,y,<j)) cos 2 d(j) 
<S>f 



dy \ r y dr y 



P(r y ,y,' 



y dy)r y dr y , 



dy > rl dr y , 



(3.1) 



(3.2) 



Iz = 





' rh{r v ) 


f"\ 


Jh{r v ) 



<t>i 



p( r y,yA) sin 2 ^>dcj> 



dy \ rl dr y 




p{r y ,y,4>) dej) 



<*(! 



y dy\r y dr y 



(3.3) 



As before, these expressions become simpler in the case of a complete revolution with azimuthal symmetry, 

rh(r v ) 



p(r y ,y) dy 



r l dr y , 



(3.4) 



1 

I x = -Iy + 2tt 



h(r y ) 



P(r v ,y) y dy 



r y dr y , 



(3.5) 



and in this case Ix = Iz-, as expected from symmetry arguments. Further, assuming p = p{r y ) the 
expressions simplifies to 



2tt 



Ix — Iz = yly 



rl p{r y ) [f 2 (r y ) - h (r y )\ dr y , 



2tt 



p{r y ) f2{r y ) ~ h(r y Y 



r y dr y 



(3.6) 
(3.7) 



We can verify again that the property Ix = Iz > Iy /2 appears in the case of azimuthal symmetry, Eq. 
(|3.5|) . This property is also satisfied in the case of incomplete solids of revolution, if conditions analogous to 
Eq. (|2.9|) for the <fi angle are fulfilled. Moreover, the theorem given by Eqs. l|2.14ll2T^|l is also held by these 
formulae, giving an alternative way for the calculation. Finally, the triangular inequalities also hold. 

These formulae are especially useful in the case in which the generating functions /i (x), f 2 {x) do not 
admit inverses, since in such case we cannot find the corresponding inverse functions gi (x) , g 2 (x) to generate 
the same figure by rotating around the X— axis. This is the case in the following example 
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Figure 5: Solid of revolution created by rotating the generating functions fa (x) = 0, and fa (x) = h + 
Asin( 1 ^) around the Y—axis. The x variable lies in the interval [0,R]. From the picture it is clear that 
one of the generators do not admit an inverse. 



Example 3 Calculate the Mi's of a homogeneous solid formed by rotating the functions 

"rvKX\ 



fl(x) = 0; fa(x) = h + Asm(^ 



R J ' 



(3.8) 



around the Y — axis (see Fig. 0), where the functions are defined in the interval x 6 [0, R], and n are positive 
integers. We demand h > \A\, ifn> 1; besides, if n = 1 and \A\ > h we demand A > 0. These requirements 
assure that fa (x) > fa (x) for all x £ [0, R\. The mass of the solid, obtained from jA.Sji) reads 



M 



irR 2 p 



nwh + 2A(-l) 



n+l 



id replacing the generating functions into the Eqs. 61 \3. r l\j we get 

R 2 / nirh + 4A(6n- 2 ir- 2 - l)(-l) r 



Ki = — 



K 2 

K x = K 2 =^L 



rnrh + 2A(-l) n+1 
1 3inrh\2h 2 + 3A 2 } 



4A{-l) n+1 [9h 2 + 2A 2 



18 [rmh + 2A(-1)"+ 1 ] 

the position of the CM (obtained from Eqs. \A.f^A.l]) , and the RG's for axes that pass through the CM read 



rcM = (0,ycjvf,0) 



Vcm 



2mr/i 2 + nnA 2 + 8Ah (-1) 



n+l 
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K 2 
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nnh + 2A(-l) 



Kz c — Kz Vcm 



Observe that fa (x) does not have inverse. Hence, we cannot generate the same figure by constructing an 
equivalent function to be rotated around the X—axis*. 

Finally, it worths pointing out that by considering homogeneous solids of revolution, we obtain the same 
expressions derived with a different approach in Ref. [Hj- 



4 Moments of inertia based on the contour plots of some figures 

Suppose that we know the contour plots of certain solid in the XY plane, i.e. the surfaces shaped by the 
intersection between planes parallel to the XY plane and the solid (see Fig. HJJ). Assume that for a certain 
value of the z coordinate, the surface defined by the contour is bounded by the functions fa (x, z) and fa (x, z) 
in the y coordinate, and by xq (z), Xf (z) in the x coordinate, as shown in the frame on the upper right 
corner of Fig. |SJ Let us form a thin plate of thick dz with the surface described above. In turn we can divide 
such thin plate into small rectangular boxes with surface dx dy and depth dz as shown in Fig. it is well 
known from the literature PP that the MI with respect to the X—axis in cartesian coordinates reads 

Ix= (y 2 + z 2 ) dm 
Jv 

'Indeed, we can find the MI of this solid by rotating around the X—axis. We achieve it by splitting up the figure in several 
pieces in the y coordinate, such that each interval in y defines a function. However, it implies to introduce more than two 
generating functions and the number of such generators increases with n, making the calculation more complex. 
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Figure 6: Contour plots of a solid utilized to calculate its moments of inertia. 

now, since our infinitesimal elements of volume are rectangular boxes with volume dV = dx dy dz, the 
contribution of each rectangular box to the MI around the X— axis is given by 

dlx = {y 2 + z 2 ) dm = (y 2 + z 2 ) p (x, y, z) dx dy dz , 

integrating over all variables, we obtain 



r*i ( r x f( z ) 


rh(x,z) 




/ V 2 P dy 


i z yj x a (z) 
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The procedure for ly and Iz is analogous, the results are. 
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Jx (z) 



fi(x,z) 

h{x,z) 

h(x,z) 



P dy 
P dy 



dz . (4.1) 

da; | z 2 dz , (4.2) 
x 2 dx \ dz . (4.3) 



Once again, we can check that results (14.11 14.21 14.3(1 satisfy Eq. (|2.16|) . This equation gives us another 
way to calculate the three Mi's. Finally, the formulae fulfill the triangular inequalities. 

An interesting case arises when we consider the Mi's of thin plates. Suppose a thin plate lying on the XY 
plane, a (x, y) denotes its surface density. This solid is generated by contour plots with volumetric density 

p(x,y,z) = cr(x,y) S (z) , (4.4) 

where 8 (z) denotes Dirac's delta function. Replacing Eq. I|4.4II into the general formula l|4.1|) we get 



Ix = 



Zo 

Xf(z) 

X (z) 



#i (z) 6 (z) dz+ H 2 (z) 6 (z) z z dz 



h{x,z) 



y 2 a(x,y) dy 



fl{x,z) 



dx ; H 2 (z) 



/(*) 

x (z) 



h{x,z) 



o-(x 7 y) dy 



fi(x,z) 



dx 



and using the properties of S (z) we get 

Ix = Hi (0) = 



f(0) 
x (0) 



h{x,0) 



y 2 a(x,y) dy 



dx 
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the z coordinate is evaluated at zero all the time, hence there is only one contour, we write it simply as 



Ix 



y <?{x,y) dy 



dx . 



Similarly Iy, Iz can be evaluated replacing Ij4.4(l into (|4.2(l and (|4.3|) 



AC*) 



a; dx ; Iz 



lY ■ 



(4.5) 



(4.6) 



Hence, Eqs. (|4.51 14. 6| give us the MPs for a thin plate delimited by f\ (x) and f 2 (a;) and by xo, x/; with 
surface density a (x,y). It worths noting that the second of Eqs. (|4.6|) arises from the application of l|4.4|l 
into (|4.3|l without assuming the perpendicular axes theorem; showing the consistency of our results^. 

The formulae shown in this section are written in terms of generating functions as in the previous 
sections. However, these generators are functions of several variables. In developing these formulae we have 
not used any particular symmetry; nevertheless, explicit use of some symmetries could simplify many specific 
calculations as shown in the following examples. 




Figure 7: Frustum of a right elliptical cone. The shadowed surfaces show a contour for a certain value of 
the z coordinate, as well as its projection onto the XY — plane. 

Example 4 Let us assume a truncated right elliptical cone as shown in Fig. ^ Such figure is characterized 
by the semi-major and semi-minor axes in the base (denoted by a\,b\ respectively), its height H, and its 
semi-major and semi-minor axes in the top (denoted by 02, b 2 respectively). Suppose that the truncated cone 
is located such that the major base lies on the XY plane and the center of such major base is on the origin 
of coordinates, as shown in Fig. ^ Now, since we are assuming that the figure is not oblique, then all the 
contours (see right top on Fig. ^ are concentric ellipses centered at the origin, with the same eccentricity. 
Therefore, it is more convenient to describe such ellipses in terms of their eccentricity e and the semi-major 
axis a (z). The contours are then delimited by 



fi{x,z) = -J a(z) -x 2 (1-e 2 ); f 2 (x,z) = J a(z) 



(!-£=») , 



x (z) 



-a (z) ; Xf (z) = a (z) 



a(z) 



(4.7) 



where /1 2 (x, z) are the functions that generate the complete ellipse of semi-major axis a (z) and eccentricity 
e (independent of z). By simple geometric arguments, we could see that the semi-major axis of one contour 
of the truncated cone at certain height z is given by^ 



a (z) = ai + 



a 2 - ai 



(4.8) 



tFor students not accustomed to the Dirac's delta function and its properties, we basically pass from the volume differential 
p dV to the surface differential a dA. 

*The semi-minor axes follow a similar equation replacing ai ; 2 — * 61, 2- From such equations we can check that the quotient 



is constant if we impose ^3- 



So the latter condition guarantees that the eccentricity remains constant. 
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Assuming that the density is constant Eq. becomes 



Iz 



'zll 

3 



H 



a(z) 
-a(z) 



\fa{z) 



2 2 
— X 1 



dx > dz 



-2pVT 



H 



(I 



C Mz) 










y a ( z Y 


— X 2 


[J-a(z) 







x 2 dx > dz 



where we have already made the integration in y. Integration in x yields 



a (z) dz 



id taking into account the Eq. we find 



Iz = (2 - e 2 ) \/l - e 2 [af + a* + a x a\ + a\a 2 + ajaj] 

Now, the mass of the figure is obtained from Eq. \A.lty and reads 

irpH 



M 



e 2 (a\ + a\ + aia 2 ) 



(4.9) 



(4.10) 



(4.11) 



therefore, the radius of gyration K\ could be written as 



Ki = — 2 
z 20 v 



af + a\ + aio% + a\a 2 + a\a\ 



aia2 + a\ + a\ 



Further, K\ and Ky can be derived from Eqs. \4-l\ \4-*\ \4-H\ ) obtaining 

3 [(l - e 2 ) (af + a\ + aia\ + a\a 2 + ajaf) + AH 2 [a\ 



K 2 



Ki = 



in 2 
6 a l 



20 



(a 2 + a\ + aia 2 ) 



-,aia 2 ) 



3 [(af + a\ + ai al + a\a 2 + a\a 2 ) + AH 2 (a 2 + \a\ + \a ia2 )] 



20 



(aj + a 2 + aia 2 ) 



when e = we get the radii of gyration of a truncated cone with circular cross section. In addition, when 
e = and a 2 = 0, the RG's reduce to the expressions for a cone with the axes X and Y passing through its 
base. Setting s = 0, a\ = 0, we also get the RG's of a cone but with the X, Y axes passing through its vertex. 
Using e — 0, a\ = a 2 we obtain the RG's of a cylinder. Finally, with a 2 = we get a cone with elliptical 
cross section, and when ai = a 2 we obtain a cylinder with elliptical cross section. Now, if we are interested 
in the MI for coordinates (Xc,Yc, Zc) passing through the CM, we should calculate the position of the CM 
from Eqs. \A . 9HA . J7)) and use Steiner's theorem obtaining 



rcM = (0,0, zcm) ; Zcm 



[2a 1 a 2 + a\ + 3a 2 ) H 
4 (aia 2 + a\ + a 2 .) 



K 



Xr 



K 2 
*^X 



z, 



cm i Ky, 



c 



K 2 



z 



CM i 



K 



-2 



K% 



(4.12) 
(4.13) 



Example 5 Frustum of a right rectangular pyramid: The contour plots are rectangles. Since the figure is 
not oblique, the ratios between the sides of the rectangle are constant. We define a\ , b\ the length and width 
of the major base; a 2 ,b 2 the dimensions of the minor base, and H the height of the solid, from which we 
have 



h 
"l 



a 2 



Hz) 
a (z) 



for all z e [0, H] , 
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where it was assumed that the major base of the figure lies on the XY plane centered in the origin with the 
lengths a\ parallel to the X—axis and the widths b\ parallel to the Y—axis. The contours are delimited by 

c c 
h(x,z) = ! h(x,z) = -a(z) , 

a(z) a(z) 
xo (z) = ^ ; x f W = ~Y~ ' 

The functional dependence on z is equal to the one in example^ so a (z) is also given by Eq. l4-8{ ). The 
integration of Eq. \4-ty gives 



I7. 



cp 
12 



11 



(l + c 2 ) / a (z) dz , 



which is very similar to Iz in Eq. \4-^ for the truncated cone with elliptical cross section, and since a (z) 
in this example is also given by Eq. \4-.3j), the result of Iz for the truncated pyramid is straightforward by 



an 



alogy with Eq. j4.1L\ ) 



Iz = 



cpH 



1 + c J (aj + a 2 + aia 2 + a*a2 + a 1 a 2 ) 



the mass is obtained from Eq. or by analogy with Eq. \4-H\j 



cpH r 2 2 n 
M = —— [a-^ + a 2 + aia 2 J 



and the radius of gyration becomes 



A'l 



20 



a\ + a\ + a\a\ + a\a2 + a 2 a\ 
a\ + a\+ a\a2 



When a,2 — we get the RG of a pyramid, if ai = a 2 we obtain the RG of the rectangular box. The RG's 
K^r, Ky are given by 

2 



K 2 



Ky 



3 T2 (£rj ( a i + a i + a i<4 + a i a 2 + ajaf) + H 2 (|aia 2 + \a\ + a|) 
5 [a\ + a\ + aia 2 ] 

3 Y2 (af + a\ + a x a\ + a\a 2 + a\al) + H 2 (±aia 2 + \a\ 
5 



2 1 a 2 ) 



[a 2 + a\ + aia 2 ] 



Finally, the expression for the position of the CM coincides with the results in example Eq. \4-l'^i) with 
the corresponding meaning of ai,a 2 in each case. The similarity of all these results with the ones in example 
^| comes from the equality in the modulation function of the contours a (z). More about it later. 



Example 6 The general ellipsoid, centered at the origin of coordinates is described by 



2 9 2 

x y z 



(4.14) 



we shall assume that a > b > c. A more suitable way to write Eq. \4-14V * s the following 

y 2 = (a{zf-x 2 ){l-E 2 ) , 

a (z) 



1- 



b(z) 
a (z) 



= \il- I - 

a 



(4.15) 
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For fixed values of z, we get ellipses whose projections onto the XY— plane are centered at the origin with 
semi-major axis a{z) and semi-minor axis b(z). The Eqs. \4-15\j show that such ellipses have constant 
eccentricity, and so we arrive to the delimited functions of Eqs. j^.7| ) with a (z) and b (z) given by Eqs. 

Therefore, the first two integrations are performed in the same way as in the truncated elliptical cone 
explained in example ^[ Then we can use the result in Eq. ( except for the limits of integration in Z ), 

the last integral is carried out by using Eqs. \4-.15\j . 



a (z) dz 



—na cp (2 



2 )v/T 



The mass of the ellipsoid is given by M = (4/3)irpabc = (4/3) itpa 2 c^f\ — e 2 , so that 



a 2 2 - e 2 



5 5 

it could be seen that the RG is independent of c, this dependence has been absorbed into the mass. Similarly, 
we can get the RG's K\ and K Y applying Eqs. \4-l\ \4- 6 ^ , the results are 



K 2 = \(b 



K 2 =\{a 2 + c 2 ) 



in this case all axes pass through the center of mass of the object. 

Observe that the Mi's for the general ellipsoid (example |5J) were easily calculated by resorting to the 
results obtained for the truncated cone with elliptical cross section (example 0}; it was because both figures 
have the same type of contours (ellipses) though in each case such contours are modulated (scaled with the 
z coordinate) in different ways. This symmetry between the profiles of both contours permitted to make 
the first two integrals in the same way for both figures, shortening the calculation of the MI for the general 
ellipsoid considerably. 

As for the truncated cone with elliptical cross section (example and the truncated rectangular pyramid 
(example EJl , they show the opposite case, i. e. they have different contours but the modulation is of the 
same type, this symmetry also facilitates the calculation of the MI of the truncated pyramid. We emphasize 
that this kind of symmetries in either the contours or modulations, can be exploited for a great variety of 
figures to simplify the calculation of their MPs. 



5 Applications utilizing the calculus of variations 



Y 

















x 


X 


) X 

f 



Figure 8: Optimization of the generating function to minimize the moment of inertia of a solid of revolution, 
the mass is the constraint and the solid lies in the interval [xQ,Xf] of length L. 

In all the equations shown in this paper, the MPs can be seen as functionals of some generating functions. 
For simplicity, we take a homogeneous solid of complete revolution around the X— axis with f\ (x) = 0. The 
MPs are functionals of the remaining generating function, from Eqs. 1(2.41 12.7(1 and relabeling /2 (x) = f (x), 
we get 

Ix[f] = %Hf(x')*dx', (5.1) 
h'lf] = I -^l+np H x' 2 f(x') 2 dx'. (5.2) 
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Then, we can use the methods of the calculus of variations (CV) *, in order to optimize the MI. To figure 
out possible applications, imagine that we should build up a figure such that under certain restrictions (that 
depend on the details of the design) we require a minimum of energy to set the solid at certain angular 
velocity starting from rest. Thus, the optimal design requires the moment of inertia around the axis of 
rotation to be a minimum. 

As an specific example, suppose that we have a certain amount of material and we wish to make up a 
solid of revolution of a fixed length with it, such that its MI around a certain axis becomes a minimum. To 
do it, let us consider a fixed interval [xo,a;/] of length L, to generate a solid of revolution of mass M and 
constant density p (see Fig. [SJ. Let us find the function / (x), such that Ix or ly become a minimum. Since 
the mass is kept constant, we use it as the fundamental constraint 

f X f 2 

M — Tip / (a;') dx' = constant. (5.3) 

J XQ 



In order to minimize Ix we should minimize the functional 

Gx [f] = f f g (/, x) dx' = I x [f] - Xnp p f (x 1 ) 2 dx' (5.4) 

J XQ J XQ 

where A is the Lagrange's multiplicator associated with the constraint (|5.3|) . In order to minimize Gx [/], 
we should use the Eulcr-Lagrange equation |5j 

SGxjf] _ dg(f,x) d dg(f,x) 

5f(x) df dxd(dffdx) 1 ' ' 

obtaining 

2npf (xf - 2nXpf (x) = , (5.6) 



Sf(x) 

whose non-trivial solution is given by 

f(x) = Vx = R. (5.7) 



Analizing the second variational derivative we realize that this condition corresponds to a minimum. 
Hence, Ix becomes minimum under the assumptions above for a cylinder of radius y/X, such radius can be 
obtained from the condition i|5.3[) . yielding R 2 = M/npL and Ix becomes 

1 M 2 

Ix, cylinder = ~ 7" ■ (^'^) 

2 npL 



Now, we look for a function that minimizes the MI of the solid of revolution around an axis perpendicular 
to the axis of symmetry. From Eqs. I|5.21 15. we see that the functional to minimize is 

Gy[f] = ^ + irp p x' 2 f {xf dx' - Xnp H f {xf dx' , (5.9) 

JXQ J XQ 

making the variation of Gy [/] with respect to / (x) we get 

/ {xf = 2 (A - x 2 ) ee R 2 - 2x 2 , (5.10) 
where we have written 2A = R 2 . By taking xq = —L/2, Xf = L/2, the function obtained is an ellipse 
centered at the origin with semimajor axis R along the Y— axis, semiminor axis R/V2 along the X— axis, 
and with eccentricity e = 1/V2- When it is revolved we get an ellipsoid of revolution (spheroid); such 
spheroid is the solid of revolution that minimizes the MI with respect to an axis perpendicular to the axis 
of revolution. From the condition l|5.3|l we find 



§The reader not familiarized with the methods of the CV, could skip this section without sacrifying the under- 
standing of the rest of the content. Interested readers can look up in the extensive bibliography concerning this topic, 
e.g. Ref. 0. 
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Figure 9: (a) Elliptical function that generates the solids of revolution that minimize ly . The shaded region 
is the one that generates the solid, (b) Truncated spheroid obtained when the shaded region is revolved. 

In the most general case, the spheroid generated this way is truncated, as it is shown in Fig. and the 
condition R > L/y/2 should be fulfilled for / (x) to be real. The spheroid is complete when R = L/y/2, and 
the mass obtained in this case is the minimum one for the spheroid to fill up the interval [— L/2, L/2], this 



minimum mass is given by 
from (O, IpTTrH) . fCT|) and (pTT!>|) we find 



^ A 1 1 i n 



npL 



*-Y, spheroid 



M min L 2 [5/i + 5u 2 - 1] 
60 



M 



M n 



(5.12) 



(5.13) 



Assuming that the densities and masses of the spheroid and the cylinder coincide, we estimate the quotients 

iY.cylinder _ {^t 1 + ~ l) 



*Y, spheroid 5/x (/i + 1) 

inder ( * J- — 2 

7 = v ^ ^ 

■* X , spheroid \ t) 



< 1, 



< l 



(5.14) 



Eqs. (|5.14|) show that Iy,sph < W,cyi while Ix.cyi < Ix,sph- In both cases if M >> M m i n the Mi's of the 
spheroid and the cylinder coincide, it is because the truncated spheroid approaches the form of a cylinder 
when the amount of mass to be distributed in the interval of length L is increased. 

On the other hand, in many applications what really matters are the Mi's around axes passing through 
the CM. In the case of homogeneous solids of revolution the axis that generates the solid passes through 
the CM, but this is not necessarily the case for an axis perpendicular to the former. If we are interested in 
minimizing Iy c , i.e. the MI with respect to an axis parallel to Y and passing through the CM, we should 
write the expression for Iy c by using the parallel axis theorem and by combining Eqs. (|5.21 15.31 IA.2|) 

Ix[f] 



Irdf] 



JY/:> 



j x' 2 f (x') 2 dx 

J Xn 



dx' 



x' f(x') 2 dx' 



thus, the functional to be minimized is 

G Yc [f]=lY c [f]-\np 



/0" 



dx' 



(5.15) 



(5.16) 



(5.17) 



after some algebra, we arrive to the following minimizing function 

/ (x) 2 = R 2 — 2 (x — xcm) 2 , 

where we have written 2 A = R 2 . It corresponds to a spheroid (truncated, in general) centered at the point 
{xcm, 0, 0) as expected, showing the consistency of the method. 

Finally, it worths remarking that the techniques of the CV shown here can be extrapolated to more 
complex situations, as long as we are able to see the Mi's as functionals of certain generating functions. The 
situations shown here are simple for pedagogical reasons, but they open a window for other applications with 
other constraints", for which the minimization cannot be done by intuition. For instance, if our constraint 
consists of keeping the surface constant, the solutions are not easy to guess, and should be developed by 
variational methods. 



11 Another possible strategy consists of parameterizing the function f (x), and find the optimal values for the 
parameters. 
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6 Analysis and conclusions 

Most textbooks report Mi's for only a few number of simple figures. By contrast, the examples illustrated 
in this work have been chosen to be more general, and can also cover many particular cases. On the other 
hand, in the specific case of solids of revolution, only the MI with respect to the symmetry axis is usually 
reported. Perhaps the most advantageous feature of the methods developed in this paper is that the three 
moments of inertia Ix, Iy , Iz can be calculated by applying the same limits of integration, and we do not 
have to worry about the partitions. It is because all three moments of inertia are written in terms of the 
generating functions of the solid. For instance, any solid of revolution acting as a physical pendulum provides 
an example in which the MI with respect to an axis perpendicular to the axis of symmetry is necessary, an 
specific case is example[H]for the Gaussian bell (see appendix IB. Moreover, we examine the conditions for 
the perpendicular Mi's to be degenerate, we find that this degeneracy occurs even for inhomogeneous solids 
as long as the density has an azimuthal symmetry Remarkably, even for incomplete solids of revolution with 
the azimuthal symmetry broken, such degeneracy may occur under certain conditions. 

Finally, we point out that for solids of revolution in which densities depend only on the height of the 
solid, the expressions for the MI become simple integrals, such fact makes the integration process much 
easier. Simple integrals are advantageous even in the case in which we cannot evaluate them analytically. 
Numerical methods to evaluate Mi's utilize typically the geometrical shape of the body; instead, numerical 
methods for simple integrals are usually easier to manage. 

As for the technique of contour plots, we can realize that many different figures could have the same type 
of contours though a different modulation of them, one specific example is the case of a cone with elliptical 
cross section and a general ellipsoid, in both solids the contours are ellipses but they are modulated (scaled 
with the z coordinate) in different ways. On the other hand, in some cases the contours are different but 
the modulation is of the same type, for example a cone and a pyramid has the same type of modulation 
(scaling) but their contours are totally different. In both situations we can save a lot of effort by making 
profit from the similarities. The reader can check the examples 01 |S] and EJ in order to figure out the way in 
which wc can exploit these symmetries in practical calculations. 

Furthermore, textbooks always consider homogeneous figures to estimate the MI, assumption that is not 
always realistic. As for our formulae, though they simplify considerably when we consider homogeneous 
bodies, the methods are tractable in many cases when inhomogeneous objects are considered, allowing more 
realistic results. 

Another interesting remark is that these methods can be used to calculate products of inertia in the case 
in which the complete tensor of inertia is necessary. Moreover, expressions for the CM proceed by similar 
arguments, and the appendix [S] shows some formulae for the CM of solids of revolution and solids built up 
by contour plots. In such expressions we realize that the same limits of integration defined for the calculation 
of the Mi's are used to calculate the CM. 

Finally, since all our formulae for Mi's depend on certain generating functions, we can see the MI of 
a wide variety of figures as a functional, making the Mi's suitable to utilize methods of the calculus of 
variations. In particular, minimization of the MI under certain restrictions is possible utilizing variational 
methods, it could be very useful in applied physics and engineering. 

The authors acknowledge to Dr. Hector Miinera for revising the manuscript. 

A Calculation of centers of mass 

A.l CM of solids of revolution generated around the X — axis 
Taking into account the definition of the center of mass for continuous systems 

/ f dm 

TCM = —f-. > 

J dm 

we can get general formulae to calculate the CM of a solid by using a similar procedure to the one followed 
to get Mi's. First of all we calculate the total mass based on the density and the geometrical shape of the 
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body. In the case of solids of revolution around the X— axis, the total mass of the solid is given by 



M 



and the CM coordinates read 



r x f ( rh(?) 








/ p(x,r x ,6) dO 


r x 


Ix \Jfi(x) 


Je 





X. 



1 f xf \ f^ 2 ^ 



CM = 



M 



x o \Jfl(x) 



p(x,r x ,9) dO 



r x dr x > dx , 



dr T > x dx 



Ycm 



^CM 



1 f Xf { ff 2 ^ 

M J X0 \J h ( x ) 



1 f f 



M 



x f \ rh( x ) 



p (x, r x ,9) cos 9 d9 



p (x, r x , 9) sin# d9 



r\ dr x > dx 



dr x > dx , 



(A.l) 

(A.2) 
(A.3) 
(A.4) 



the limits of integrations are the ones defined in Sec. 12.11 In the case of complete revolution with p = p(x,r x ), 
we obtain Ycm = Zcm = 0, as expected from the cylindrical symmetry. 

A.2 CM of solids of revolution generated around the Y — axis 

By using the coordinate system and the limits of integration defined in Sec. 01 we can evaluate the CM for 
solids of revolution generated around the Y— axis obtaining 



X, 



1 f Xf \ [f 2 ( ry ^ 



CM 



M 



Y 



1 f Xf { f^ 2 ^ 



CM 



^CM 



M 



M 



x o \ J h(r y ) 



1 f xf { f f2 ^ Ty ^ 

M Jx \Jh(r y ) 



<$>} 
<t>o 

4>f 

<t>0 

<t>f 

<t>0 



$m(4>)p{r y ,y,4>) d(j) 



dy > r y dr y 



p(r y ,y, 



y dy\r y dr y 



p(r 3/ ,y,0)cos(0) d(j> 



dy } r v dr. 



y <*iy , 



x a \Jh{r y ) 



4>f 



00 



p(r v ,y,<f>) dej) 



dy } r y dr y 



(A.5) 
(A.6) 
(A.7) 
(A.8) 



In the case of a complete revolution with p = p(r y ,y), we get Xqm = Zcm = 0, due to the cylindrical 
symmetry. 

A.3 CM of solids formed by contour plots 

In this case, we use the limits of integration defined in Sec. 0] The CM coordinates read 



I r z f ( r x f( z ) 



XCM 



ycM 



M 



z o I Jxo(z) 



I r z i I r x f( z ) 

MJz \Jx (z) 

I r z f ( r x f( z ) 



ZCM = M 



M 



Z I J Xq[Z 



h{x,z) 

fl( x ,z) 
fl{x,z) 

fl(x,z) 



p(x,y,z) dy 



x dx > dz , 



p(x,y,z) y dy 



dx } dz 



h( x ,z) 





rx f (z) 


rf2(x,z) 




_ Jx (z) 


Jh{x.z) 



p(x,y,z) dy 
p(x,y,z) dy 



dx > z dz 



dx > dz 



(A.9) 
(A.10) 
(All) 
(A.12) 
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B Some additional examples for the calculation of moments of 
inertia 

In this appendix we carry out additional calculations of moments of inertia for some specific figures, by 
applying the formulae written in sections [3 [21 and in order to illustrate the power of the methods. 

B.l Examples of moments of inertia for solids of revolution 




Figure 10: On left: cylindrical wedge generated around the Y—axis. On right: A bell modelated by two Gauss' 
distributions rotating around the Y— axis, the bell tolls around an axis perpendicular to the sheet that passes 
through the point B 

Example 7 MI for a cylindrical wedge (see Fig. \1U\) . Let us consider a cylinder of height h and radius b, 
whose density is given by 

p, a < <b < 2ti — a 



o. 



a < <t> < a 



and that is generated around the Y—axis by means of the functions f\ (x) = and f 2 (x) = h. From Eq. 
ifAgj) and Eqs. yj.ll'J.tf) we find 

M = (n - a)phb 2 , 



M 



b 2 f sin 2a 
~8 V (tt - a) 



3 



h 



Mb 2 



I9. 



M 



b 2 ( sin 2a 
~8 \ + (tt - a) 



3 



In order to calculate the moments of inertia from axes passing through the center of mass we use Eqs. 
li.Mjl) to get 

h 2 b 

X C m = ; Y CM = - ; Z CM = sin a. 

2 3 (it — a) 

Example 8 Mi's for a Gaussian Bell. Let us consider a homogeneous hollow bell, which can be reasonably 
described by a couple of Gaussian distributions (see Fig. \10\) . 

h{x) = Ae- a * 2 ■ f 2 {x)^Be-^ 2 , 
where a,[5,A,B are positive numbers, A < B, and a > (3. The Mi's are obtained from \3.(k 



ly = irp 



B A 

~f¥~a^ 



Ix=Iz = -^ 



up 


B A' 


irp 


B 3 


A 3 ' 


T 


J 1 ~ ^_ 


+ i 


. 


a 



Besides, the mass and the center of mass position read 



M = up 



B_A 

J~a 



Ycm = T 



aB 2 - (3A 2 
aB - j3 A 



Xcm — Zcm = 
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When the bell tolls, it rotates around an axis perpendicular to the axis of symmetry that passes the top of the 
bell. Thus, this is a real situation in which the perpendicular MI is required. On the other hand, owing to 
the cylindrical symmetry, we can calculate this moment of inertia by taking any axis parallel to the X—axis. 
In our case the top of the bell corresponds to y = B, and using Steiner's theorem it can be shown that 



I 



I x + MB(B - 2Y C 



CM) 



Ix.B = 



71 P 



18a 2 /3< 



[a 2 B(9 + llB 2 f3) + f3 2 A{9ABa - 2A 2 a - lSB 2 a - 9)] 



B.2 Examples of Mi's by the method of contourplots 

Example 9 A thin elliptical plate: for this bidimensional object, we can use Eqs. j4-5\ \4--°M , ^ e delimited 
function can be taken from J^.7| ) but with 2 = 0. The results are 



Iy = M-, 



M- 



M = rraab . 



once again, these axes passes through the center of mass of the figure. 



Y 




X 

Figure 11: Arbitrary quadrilateral, the dimensions are indicated in the drawing. 



Example 10 An arbitrary quadrilateral, (see Fig. Ml)) : This is a bidimensional figure, so we apply Eqs. 
\4-5[ \4-$il - The bounding functions are given by 

fi {x) = 0, 

if < x < ai 

fe-M + M°i+« a )-fr.°i if ai <x< ai +a 2 

— + ^ (ai + a-2 + a 3 ) if a\ + a 2 < x < ai + a 2 + a 3 
and the moments of inertia read 



^x 

ai 



h{x) 



Ix= V2 ^ ^ + h ^ ^ + H ^ + ^ + /l ^ 3 ] 



(B.l) 



[l2aia 2 a 3 /i 2 + (Aa\a 2 + 3a 3 + a 2 ) h\ + (4a 2 a 3 + 3a 



+6a 2 a 2 [hi + h 2 ) + ^a\h 2 {2a\ + a 3 ) + 6a 3 /i 2 (a\ + a^)] 



(B.2) 



(B.3) 
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the center of mass coordinates are given by 



xcm = [3aia 2 (hi + h 2 ) + 3a 3 /i 2 (01 + a 2 ) 

+hi (2a\ + a%) + h 2 (2aj + a§)] 



(B.4) 



Vcm = \a2h1h1 + (ai + a 2 ) h\ + (a 2 + a 3 ) hi] 



(B.5) 



with 



M = — [aihi + a 2 (hi + h 2 ) + a 3 h 2 ] 



(B.6) 



The quantities Ix, He Mi an d M; are invariant under traslations in x, for example Ix might be calculated 
as 



where we have performed a traslation Ax to the left. It is equivalent to make the change of variables 
u = x — Ax. This property can be used to evaluate the integrals easier. Specifically, the piece of Fig. luina 
in the interval ai < x < ai + a 2; can be traslated to the origin by using Ax = ai; and the piece of this figure 
lying at ai + a 2 < x < ai + a 2 + ct3 can be also traslated to the origin with Ax = ai + a 2 . On the other hand, 
though the quantities Iy and xcm ore not invariant under such traslations, the same change of variables 
simplifies their calculations. This strategy is very useful in solids or surfaces that can be decomposed by 
pieces (i.e. when at least one of the generator functions is defined by pieces). For example, the same change 
of variables could be used if we are interested in the solid of revolution generated by the surface of Fig. 
Finally, from Eqs. \B. WB. 6)) we can obtain many particular cases, some of them are 



• 


hi 


= h 2 


trapezoid. 


• 


ai 


= hi 


= arbitrary triangle. 


• 


a 2 


= 0, 


o>l = 0,3 = ^7= = ^ = i£ equilateral triangle of side length L 


• 


«i 


= hi 


= 0, a 2 = 03 = ^ = -j equilateral triangle of side length L. 


• 


«i 


= a.3 


= hi = triangle with a right angle. 


• 


h 2 


= hi 


, a 2 = arbitrary triangle. 


• 


ai 


= a 3 


= 0, h2 = hi rectangle. 



C Table of moments of inertia 

In tabic □ on page 1221 the moments of inertia for a variety of solids of revolution generated around the 
Y~ axis are displayed, such table includes the function generators and any other information necessary to 
carry out the calculations by means of our methods. The surfaces that generates the solids are displayed in 
Fig. ^] Observe that the first of these surfaces generates a torus with elliptical cross section and the second 
one generates a truncated hollow cone. 

Finally, there are some conditions for certain parameters of these figures. In Fig. (d), a > 0, and b > 2a; 
in Fig. (e) b > 0, and a > 0, in Fig. (f), n > 0; for Fig. (g), n > —2/3; in Figs, (c) and (h) a can also be 
negative. 
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c) 





Figure 12: Surfaces that generates the solids whose moments of inertia appears on the tabled 
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